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ABSTRACT 

We use the BzK color selection proposed by iDaddi et al. ( 2004 ) to obtain a sample of 1092 
faint star-forming galaxies (hereafter sBzKs) from 180 arcmin 2 in the Subaru Deep Field. This 
sample represents star-forming galaxies at 1.4 < z < 2.5 down to Kab — 23.2, which roughly 
corresponds to a stellar-mass limit of ~ 1 x 10 10 M Q . We measure the angular correlation function 
(ACF) of these sBzKs to be w(6) = (0.58 ± 0.13) x 6»["]-°- 8 and translate the amplitude into the 
correlation length assuming a reasonable redshift distribution. The resulting value, r Q = 3.2+"'^ 
hr 1 Mpc, suggests that our sBzKs reside in haloes with a typ ical mass of 2.8x10 11 M ( t ) . Combining 
this halo mass estimate with those for brighter samples of Kong et al. ( 20061 ). we find that the 
mass of dark haloes largely increases with K brightness, a measure of the stellar mass; the dark 
halo mass increases by as much as 10 2-3 as K brightness increases by only a factor of ~ 10. 
We also find that the halo occupation number, the number of galaxies hosted in a dark halo, is 
higher for brighter sBzKs. Comparison with other galaxy populations suggests that faint sBzKs 
(Kab < 23.2) and Lyman Break Galaxies at z ~ 2 are similar populations hosted by relatively 
low-mass haloes, while bright sBzKs (Kab < 21) reside in haloes comparable to or more massive 
than those of Distant Red Galaxies and Extremely Red Objects. Using the extended Press- 
Schechter formalism, we predict that present-day descendants of haloes hosting sBzKs span a 
wide mass range depending on K brightness, from lower than that of the Milky Way up to those 
of richest clusters. 



Subject headings: cosmology: observations 
high-redshift — galaxies: photometry 



galaxies: evolution — galaxies: formation — galaxies: 



1. Introduction 

Galaxy mass is probably the most important 
parameter among those governing the evolution 
of galaxies. Observations of the local universe 
show that more massive galaxies tend to have 
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older stellar populations, lower star-formation ac- 
tivities, earlier m orphological types, a nd higher 
metallicities (e.g.. IGallazzi et"aL 12005 ). In Cold 
Dark Matter (CDM) universes, dark haloes grow 
with time through merging of less massive haloes, 
implying that haloes with different masses have 
different evolutionary histories. Observing galax- 
ies over a wide mass range is thus crucial to place 
strong constraints on galaxy evolution. 

It is extremely difficult to measure dark-halo 
masses of high-redshift galaxies. One can, how- 
ever, infer them from galaxy clustering on the as- 
sumption that galaxies reside in dark haloes, since 
CDM models predict a monotonic correlation that 
more massive haloes are clustered more strongly 
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( Mo fc White! 2002 ) . Measuring galaxy clustering 
requires a large sample from a wide area. 

Recent observations suggest that the era of z ~ 
2 is important in galaxy evolution for various rea- 
sons; the cosmic star formation rate begins to drop 
at z ~ 1 — 2 from a flat p lateau at higher red shifts 



(|Dickinson et all 120031 : iFontana et al.l 12003 ): the 



morphological type mix of field galaxies changes 
rema rkably at z ~ 1 - 2 (IKaiisawa fc Yamada 



200 if ): the number density o f QSOs has a peak at 



2 ( Richards et al. 20061 ). Various color selec- 



tion methods such as the L yman-break technique 
( Guhathakurta et al. 1990h have been proposed so 
far, in order to construct large samples of galax- 
ies at high redshifts from imaging data. However, 
color selections for galaxies at z ~ 2, the 'red- 
shift desert', have not been as successful as those 
for other redshifts, mainly because of the lack of 
distinctive spectral features in the optical wave- 
length. 

Recently, a two-color select ion with B — z' and 
z' — K has been proposed by iDaddi et all (|2004 ) 
to select galaxies at 1.4 < z < 2.5 effectively with 
a low contamination. Galaxies selected with this 
method, referred to as BzKs, make up a good ap- 
proximation of a sample of all galaxies at z ~ 2 
down to a limiting magnitude in the K band, 
roughly equal to a stellar mass limit. Two types 
of BzKs, sBzKs and pBzKs, are defined accord- 
i ng to the position in the B — z' vs z' — K space 
(Daddietal ]|2004h : V and 'p' denote star-forming 
and passive, respectively. It has been shown that 
most star-forming galaxies at 1.4 < z < 2.5 meet 
the criteria for sB zKs irrespective of the amount 



of dust extinction (jDaddi et al.l 120041 ) 



Previous observations of BzKs have been lim- 
ited primarily to bright objects with K < 22. 
Those studies have shown that bright sBzKs have 
star formation rates (SFRs) of ~ 10 2 M© yr -1 , 
stellar masses of 1Q 11 M GU and dus t extinctions 
of E(B - V) ~ 0.4 (IDaddi et all I2004L 120051: 



Kong et al.l I2006L hereafter K06). Bright pBzKs 
are found to have stellar masses of > 1Q u Mq with 
ages of ~ 1 Gyr. Clustering has also been exam- 
ined recently using samples of bright BzKs from 
wide- field surveys (K06); both pBzKs and sBzKs 
are found to be clustered strongly. From these 
studies, bright BzKs are inferred to be progeni- 
tors of the present-day massive galaxies. However, 
properties of faint, or low-mass, BzKs remain un- 



revealed due to the difficulty in obtaining deep and 
wide-field NIR images. 

Steidel et al.l (|2004f ) have proposed another 
method to select gala xies at z ~ 2 (see also 



Adelberger et al.l 12004) . This method uses U,G 



and R photometry to select far-UV bright star- 
forming galaxies, referred to as BXs/BMs. Clus- 
tering properties of BXs/BMs have been s tudied 
using large samples ( Adelberger et al.ll2005l ). since 
NIR imaging is not required. Note, however, that 
this method will miss dusty galaxies and passive 
galaxies, both of which are faint in far-UV wave- 
lengths. 

We have been conductin g a deep survey of BzKs 
in the Subaru Deep Field (jKashikawa et al.ll2004 : 



Maihara et al.ll200lh to study properties of faint, 



or low-mass, BzKs. Details of our survey and the 
photometric properties of BzKs detected from our 
survey, such as number densities, star- formation 
rates, and stellar masses, will be presented else- 
where (Motohara et al. in prep). 

In this paper, we report results of clustering 
analysis of sBzKs detected in our initial data set 
over 180 arcmin 2 . Our study provides the first 
measurements of angular correlation of faint [K < 
23.2) sBzKs. Unfortunately, the number of pBzKs 
in the present sample is not large enough to obtain 
reliable measurements of their angular correlation. 
Combining our result with those for bright [K < 
22) sBzKs given in K06, we examine properties of 
dark haloes hosting sBzKs over a wide range of K 
magnitude, or, equivalently, stellar mass. Dark- 
halo masses are inferred from observed correlation 
lengths through an analytic model for the spatia l 
clustering of dark haloes fe.g. jMo fc Whitej|2002l) . 
We also compare the dark-halo masses of sBzKs 
with those of other populations such as Lyman- 
break galaxies (LBGs). 

The structure of this paper is as follows. The 
optical and near-infrared data used to select BzKs 
are described in In [J21 we construct a sam- 
ple of faint sBzKs and check its consistency by 
comparing it with other samples. We measure the 
angular correlation function of sBzK in our sam- 
ple in []4l Then, for both our sample and K06's, 
the amplitude of the angular correlation function 
is transformed into the spatial correlation length. 
Results and discussion are presented in Ej5j and 
conclusions are given in S|6] 
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Throughout this paper, magnitudes are in the 
AB system, and we adopt cosmological parame- 
ters of f2 m o = 0.3, Qao — 0.7. and erg = 0.9. 
We assume h — 0.7 to compute the power spec- 
trum of the cosmological density fluctuations and 
physical quantities except the correlation length, 
where h is the Hubble constant in units of f 00 km 
s _1 Mpc . The correlation length is expressed in 
units of /i -1 Mpc, since most of the previous stud- 
ies adopted h = 1. 

2. Data 

2.1. Optical Data 

The Subaru Deep Field (SDF) is a blank field 
centered on (13 h 24 m 38!9, +27°29'25'.'9) (J2000). 
The SDF has deep and wide- field Subaru/Suprimc- 
Cam data of seven optical bandpasses, B, V, R, i', 
NB816, and NB9 21, obtained for the Su baru 



Deep Field Project (jKashikawa et al. 2004). In 



this paper, we combine the B and z' data with the 
K data described below, to select BzK galaxies. 
The exposure time and the 3<7 limiting magnitude 
on a 2" aperture are 595 min and 28.45 mag for 
B and 504 min and 26.62 mag for z', respectively. 
Both images have been convolved to a seeing size 
of I'.'U (FWHM), the value for the final J and 
K images (see below), and have an identical sky 
coverage of 29'7 x 36f 7 = 1090 arcmin 2 with a 
pixel scale of 0'.'202 pixel" 1 . 

2.2. Near-Infrared Data 

We observed the SDF in the J and K bands 
with the wide-field camera (WFCAM) on the 
United Kingdom Infra- Red Telescope (UKIRT) on 
2005 April 14 - 15. WFCAM is composed of four 
2048 x 2048-pixel detectors with a large spacing of 
12.3 arcmin. Each detector covers ~ 13' x 13' of 
sky with a pixel scale of C'4. We made only a sin- 
gle pointing for each passband to go significantly 
deeper than the previous wide-field K data used 
to study BzK galaxies. 

Figure [T] shows the sky area of the SDF imaged 
with Suprime-Cam and WFCAM. The Suprime- 
Cam 297 x 367 field is outlined by the thin 
solid line. The four dotted rectangles correspond 
to the regions imaged with the four detectors of 
WFCAM. The large spacing between the detec- 
tors prevented us from covering the Suprime-Cam 
field effectively. The total area overlapping with 
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Fig. 1. — The sky area of the SDF observed with 
Suprime-Cam and WFCAM. The thin solid lines 
outline the 29'7 x 36'7 area imaged with Suprime- 
Cam (B,z r ). The four dotted rectangles corre- 
spond to the four fields with WFCAM J and K 
data, we used in this study only the two subre- 
gions enclosed with thick solid lines in the north- 
west and southwest fields. 



the Suprime-Cam field amounts, however, to 410 
arcmin 2 , large enough to measure the angular clus- 
tering of faint BzK galaxies. 

The exposure time was 150 min for J and 294 
min for K . The variations in sensitivity and PSF 
among the four detectors were found to be small. 
The seeing size in the final images varied from 1'.'04 
to T/14 in K. 

Data reduction was made in the standard man- 
ner for NIR imaging data, but special care was 
taken to exclude spurious objects due to crosstalk 
of bright objects. All the images, except the one 
with the worst seeing size, were smoothed with a 
Gaussian kernel so that their PS F sizes be 1'.'14 
FWHM. The 2MASS catalog (|Skrutskie et al 



2006) was used to conduct astrometry and to de- 
rive the magnitude zero point. 

We measure the 5<r limiting magnitude on a 
2"-diameter aperture for each detector. The shal- 
lowest (brightest) limiting magnitude among the 
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four detectors is 23.3 mag in J and 23.5 mag in 
K. The variation in limiting magnitude among 
the four detectors is less than 0.2 mag for both 
bandpasses. Details of our observation and data 
reduction will be described elsewhere (Motohara 
et al. in prep). 

2.3. /f-Selected Catalog 

Object detection an d photometry are mad e 
using SExtractor 2.1.3 (jBertin fc Arnoutd 1 19961) . 
The X-band is chosen as the detection band. For 
each object detected in the K image, 2 "-diameter 
magnitude is measured for B and z' using double- 
image mode of SExtractor to derive colors. If an 
object is undetected (fainter than the 2a magni- 
tude) in B or z', its magnitude is replaced with the 
2a magnitude. We limit our catalog to K = 23.5, 
the shallowest 5er limiting magnitude among the 
four detectors. We have in total 8884 objects in 
the overlapping fields (see Figure [T]). The cata- 
log thus constructed is then correc ted for Galac- 
tic ext inction using the dust map of Schle gel et al.l 
(|l998l ): A{B) = 0.07 mag, A(z') = 0.02 mag, and 
A{K) = 0.01 mag. 

2.4. Detection Uniformity 

Clustering analysis requires high uniformity in 
object detection over the image for all passbands, 
since a variation in detection completeness may 
produce spurious clustering signals. While the 
WFCAM J and K images are very uniform, sky 
noise at the four edges of the Suprime-Cam images 
is systematically large and non-uniform because of 
dithered pointing with large angles. 

We examine detection uniformity for each of the 
four overlapping fields by dividing the field into 
a few dozen of sub-areas and estimating the sky 
noise in each sub-area for B, z', and K . On the 
basis of these noise maps, we find that the two 
subregions surrounded by thick solid lines in Fig- 
ure [T] have a good detection uniformity in all three 
passbands. These subre gions are well within th e 
good S/N area defined in lKashikawa et al.l (|2004l ). 
In what follows we call these two subregions the 
NW and SW subregions. The area of NW and 
SW subregions are 100 arcmin 2 and 80 arcmin 2 , 
respectively. 

In this paper, we examine angular clustering 
only in the NW and SW subregions. Among the 



8884 objects found in the four overlapping fields, 
3961 are included in these two subregions. The 
two eastern fields have also sub-areas with a good 
detection uniformity, but they are all too small for 
clustering analysis. 

3. BzK Galaxies in the Subaru Deep Field 
3.1. Selection of BzK Galaxies 

In the B — z' vs z' — K plane, sBzKs are de- 
fined to be objects with (z 1 — K) — (B — z') > 
—0.2, while the criteria for pBzKs are defined as 
') < -0.2 n z - K > 2.5 
Stars and foreground galax- 



(z 1 - K) - 


B - z 


Daddi et al. 


2004) 



ies are separated well fr om BzK galaxies in this 
plane (jDaddi et al.ll20041 h We also regard objects 
with K < 18 n CLASS_STAR > 0.9 as stars, 
where CLASS_STAR is the stellarity index param- 
eter from SExtractor. Although the validity of the 
BzK selection for objects fainter than K ~ 22 
has not been examined well, we assume that the 
selection and the classification into two types are 
applicable to our objects as well. 

In the selection of BzKs, special care is required 
for objects undetected in B or z', i.e., B > £?2<riim 
or z' > z' 2aUm , where B 2tT \im and z' 2a . Km are the 
2a limiting magnitudes. We treat these objects as 
follows: 

• Objects undetected in B but detected in z'. 
If their I^o-iim — z' and z' — K colors meet the 
pBzK criteria, they are classified as pBzKs, 
since their true B — z' colors, which are red- 
der than S2criim— z', also meet (z 1 —K) — (B— 
z') < —0.2. On the other hand, classification 
is not applied to those whose Bi a u m — z 1 and 
z' — K colors satisfy the sBzK criterion, since 
their true B — z' colors may be too red to be 
classified as sBzKs. 



• Objects detected in B but undetected in z'. 
If their B - z' 2r7lim and z' 2aUm - K colors 
meet the sBzK criterion, they are classified 
as sBzKs, since they do not go out of the 
sBzK region even when true z' magnitudes, 
which arc fainter than z' 2alim , are used in- 
stead. On the other hand, classification is 
not applied to those with B — z' 2lTVlm and 
z' 2cr y lm — K satisfying the pBzK criteria, since 
their true B — z' may be too blue to be clas- 
sified as pBzKs. 
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• Objects undetected in both B and z' . 

These objects are not classified, because 
their B — z' colors can take any values. 

The numbers of sBzKs and pBzKs selected are 
1092 and 56, respectively (Table HJ. Figure 
shows the B — z' vs z' — K distribution of all ob- 
jects in the NW and SW subregions. 

Our B, z', and K bandpasses are not exactly 
the same as those origi nally used to defin e the 



Bz K selection criteria in Daddi et al. (2004). We 



Gunn fc Strvken (|l983f) 's stellar spectrophoto- 



use 

metric atlas to find that the B — z' colors of stars 
defined in our syst em can be bluer than those in 
Daddi et al.l (|2004l ) up to 0.2 mag, depending on 
the spectral type, while there is little difference in 
z' — K. We also make a similar evaluation using 
the spe ctral templates of loca l galaxies (E to Im) 
given in Coleman et al. ( 198dh redshifted to z ~ 2. 
It is found that the B — z' colors of these galaxies 
in our system a r e at m ost 0.1 mag bluer than those 
in iDaddi et al.l (2004) while our z' — K colors are 



at most 0.1 mag redder. 

These results imply that applying the original 
sBzK boundary, (z' - K) - (B - z') > -0.2, to 
our data will select additional objects near the 
boundary which, if measured in the system of 



Daddi et al.l (|2004f ). may not be selected. Our 
sBzK sample will be reduced by ~ 15 % when the 
selection boundary is tightened by 0.2 mag (i.e., 
(z 1 -K)-{B- z') > 0.0). This reduction should 
be an upper limit to the effect of the bandpass 
differences, since the net offset in the B — z 1 vs 
z' — K plane of sBzKs due to the bandpass dif- 
ferences will be on average more modest than 0.2 
mag. Because this reduced sample gives a very 



similar ACF, we have decided not to adjust the 
selection boundary in our study. 

We find that the 2"-aperture K magnitudes of 
our BzKs are on average 0.3 mag fainter than the 
total magnitudes (we adopt MAG_AUTO magni- 
tudes for total magnitudes). For this reason, we 
regard that our BzK sample is flux-limited to a 
total magnitude of K = 23.2(= 23.5 - 0.3). 




Fig. 2. — B~z' vs z' —K distribution of all objects 
in the NW and SW subregions. For objects unde- 
tected in B or z' , B2 a \\m and z' 2alim have been 
used, respectively, in the plot. Filled circles rep- 
resent BzKs; those with a horizontal arrow are 
undetected in B. Crosses represent unclassified 
objects. 



Table 1: Numbers of sBzKs, pBzKs and unclassi- 
fied objects 



Detection* 


sBzK 


pBzK 


unclassified 


B z' 








o o 


1092 


40 





X o 





16 


20 


o X 











X X 








9 


total 


1092 


56 


29 



* o and x represent detection and non-detection, 
respectively. 



3.2. Number Counts 

Figure [3] compares the number counts of sBzKs, 
pBzKs, and all galaxies from our data with those 
given by K06. Since K06's data are from two wide 
fields, EIS Deep3A field and Daddi field, we plot 
the counts from these two fields separately. It is 
found that our number counts of sBzKs and all 
galaxies are consistent with those of K06 while 
pBzK are lower by a factor of 2 - 3. This dis- 
agreement in the number counts of pBzK may be 
partly due to cosmic variance. Indeed, even in the 
two fields surveyed by K06, both of which are sev- 
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eral times larger than our total area, the counts 
are different significantly. A small part of the dis- 
agreement may also result from our special treat- 
ment in the BzK selection for objects undetected 
in B or z' described in Section [57X1 which has not 
been considered in K06. 




method (|Ling et al.|[l98l iMo et al.lll992l) . Then 
the measurements for each subregion are corrected 
for the integral constraint, on the assumption that 
the true ACF obeys a power law of w{9) oc 6*~ 8 . 
Finally, the measurements from the two subre- 
gions are averaged to obtain the ACF of sBzKs 
in the SDF (Figure HJ). 

The amplitude of the ACF, A, is obtained from 
a fit of w(9) = A6~ - 8 to the data, and is given 
in Table [1 Our sBzKs have A = 0.58 ± 0.13 
[arcsec 8 ], where the errors correspond to the 
range in which the increase in \ 2 from the best-fit 
value is less than unity. Our ACF measurement 
for sBzKs could suffer from cosmic variance, since 
each of the two sub-fields is only about 10 Mpc a 
side at z = 2. A survey of a larger area is clearly 
required to obtain results robust against possible 
variance. However, we infer that the amplitude 
of the cosmic variance in our sample will not be 
intolerably large, since the ACFs from the two 
sub-fields agree within the statistical errors, and 
since the number counts of our sBzKs are consis- 
tent with those of K06. 



Fig. 3. — Number counts of sBzKs, pBzKs, 
and all galaxies as a function of total magnitude 
(MAG_AUTO). The circles, triangles, and squares 
represent sBzKs, pBzKs, and all galaxies, respec- 
tively. Our data are connected by the solid line, 
while K06's are connected by the dotted line. The 
K06 data are plotted separately for the two survey 
fields. 



4. Clustering Analysis 
4.1. Angular Correlation 

We derive the angular correlation function 
(ACF) separately for the sBzKs and pBzKs. We 
do not, however, show the results of the pBzKs, 
since reliable measurements are not obtained due 
to their small sample size (56 objects). 

We first measure the ACF of sBzKs for the 
NW and SW subregions independently with an- 
gular bin s of logflH = . 4, usi ng the estimator 



given in lLandv & Szalavi (|l993h . Random data 



with 50000 points are used. Errors in the measure- 
ments are estimated by the bootstrap resampling 




10 100 
0[arcsec] 



Fig. 4. — The angular correlation function of 
sBzKs. The error bars correspond to the la errors 
obtained from bootstrap resampling. The dotted 
line is the best-fit power law of A8~ 8 . 
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4.2. Bright Samples by K06 

Figure [5] shows the observed ACF amplitude 
against limiting K magnitude. The filled circle 
represents the amplitude of our sBzKs, while the 
six filled triangles correspond to the amplitudes 
for the six samples of K06 with different limiting 
magnitudes over K = 20.4 and 21.9. The ACF 
amplitude of sBzKs is found to decrease rapidly 
with decreasing K brightness. The brightest sam- 
ples of K06 have A ~ 15 - 20 arcsec ' 8 , which is 
about 30 times higher than ours. 

In the next subsection, the spatial correlation 
lengths, ro, are calculated from the ACF ampli- 
tudes to discuss the clustering of sBzKs in real 
space. For the reader's reference, the three shaded 
regions in Figure [5] indicate, respectively, the ex- 
pected ACF amplitude for r = 3, 10 and 20 
[/i _1 Mpc] calculated using the redshift distribu- 
tion for sBzKs adopted in the next subsection 
(Gaussian with z c = 1.9 and o~ z — 0.35); for each 
ro, the spread in amplitude corresponds to the 
adopted uncertainty in the o~ z of ±0.1. 

4.3. Spatial Correlation Length of sBzKs 

The correlation length, ro , for sBzKs can be cal- 
culated from the a mplitude of the ACF, A, usi ng 
Limber's equation ( Limber 1953 : Peebles! 1980l ) if 
their redshift distribution is known. 

4-3.1. Redshift Distribution of sBzKs 

The redshift distribution, N(z), of sBzKs has 
not been established even for bright objects with 
K < 22 owing to the lack of a large, flux-limited 
spectroscopic sample. In this study, we assume 
that the N(z) of sBzKs is not dependent on K 
magnitude, and adopt a Gaussian distribution as 
a simple approximation. We determine the central 
redshift, z c , and the standard deviation, o~ z , of the 
Gaussian as follows. 

First, we use a preliminary result of the N (z) of 
sBzKs (81 objects, most of which are K = 21 — 22) 
obtained from a X-selected spectroscopic survey 
in the EIS Deep3A field with VLT/VIMOS down 
to K — 22.0 (E. Daddi, private communication) 
to determine z c and infer o~ z . This is probably the 
most reliable measurement of the N(z) of sBzKs 
at present because it is based on the largest spec- 
troscopic sample. We fit a Gaussian to this N(z) 
after removing a strong spike at z ~ 1.5 (which 



o h 



# our sample 
A K06 





20 20.5 21 21.5 22 22.5 23 23.5 24 24.5 
Limiting K AB Magnitude 

Fig. 5. — Amplitude of the ACF, A, as a function 
of limiting K total magnitude. The filled circle 
shows the measurement for our sample, while the 
filled triangles correspond to the six samples with 
different limiting magnitudes over K — 20.4 and 
21.9 given in K06. The amplitudes for r = 20, 10 
and 3/i _1 Mpc calculated using the adopted red- 
shift distribution for sBzKs (Gaussian with z c = 
1.9 and a z = 0.35) are shown by the shaded re- 
gions, whose width corresponds to the adopted un- 
certainty in a z of ±0.1. See subsection 14.3.11 for 
details. 



could be due to large-scale structure) and limiting 
the redshift range to 1.0 < z < 3.0, and obtain z c 
= 1.9 and o~ z — 0.4. From this result, we adopt z c 
= 1.9 for our Gaussian distribution. Note that the 
Limber transformation is insensitive to the change 
in z c ; for example, changing z c by ±0.2 around 1.9 
changes r only by 2.5 %. 

Next, we determine a z and express the uncer- 
tainty in our Gaussian distribution by the range 
of a z . We adopt a z = 0.35 ±0.1 for the follow- 
ing reasons. We find that, for a fixed A value, 
the range of ro corresponding to the range 0.25 < 
a z < 0.45 covers not only the ro value calcu- 
lated using the raw N(z) from the EIS Deep3A 



field, but also from tho s e give n in iDaddi et al 
(l2004h and iReddv et al] (|2005fh Here, the data 
of iReddv et al. (2005) are another measurement 
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of the redshift distribution based on a large sam- 
ple, although the sample is not isT-limited but opti- 
cally selected and thus could be biased toward UV 
bright objects. We also find that the best-fit Gaus- 
si an parameters for t he re dshift distribution data 
of iDaddi et all (|2004l) and lReddv et all (|2005l ) are 
similar to the values we adopt above. Figure [6] 
compar es the adopted Gau ssian function with the 
N(z) of IDaddi et all (|2004h . 

It is true that a Gaussian distribution is not an 
excellent approximation to the existing N(z) data, 
especially for the preliminary result of Daddi et 
al. (private communication) with a large spike. 
However, it will be reasonable to expect that the 
true N(z) derived from a sample large enough to 
smooth out statistical noise and cosmic variance is 
not so far from a Gaussian centered at z ~ 2, since 
the isolation of 1.4 < z < 2.5 galaxies from others 
using the BzK selection will not be perfect. It may 
be worth noting that Lyman-break galaxies, which 
are also selected in a two-color plan e, are known to 



have a Gaussian-like distribution (jYoshida et al 
2006). Our purpose is not to infer the true dis- 



tribution but to find a reasonable expression of 
it which can be used in the Limber transforma- 
tion. In this sense, the adopted Gaussian and 
uncertainty, z c = 1.9 and a = 0.35 ± 0.1, are a 
reasonable parameterization, since they include all 
the variations in the existing redshift distribution 
measurements in terms of resultant ro values. It 
is also found that the ro value calculated using 
a top-hat redshift distribution of 1.4 < z < 2.5 
(ideal selection function for BzKs) falls within the 
ro range defined by the adopted uncertainty in the 
Gaussian. 

Finally, we should note that our assumption 
that the N(z) does not depend on K magnitude 
must be tested using a future spectro scopic follow- 
up, a lthough the observation of iReddv et al 
(|2005h may support this assumption. They de- 
rived N(z) for K < 21.8 and K > 21.8 (with a 
limiting magnitude of K = 24.3) separately, which 
are found not to be largely different. 

4-3.2. Spatial Correlation Length 

The correlation length of our sBzKs derived in 
this manner is ro = 3.2^o 7 [/i _1 Mpc] in comoving 
units (see also Table [2]). The errors in ro are cal- 
culated by summing up in quadrature the errors 
from the measurement of A and the errors due to 



6 - 



2 - 




Fig. 6. — The redshift distribution of sBzKs. The 
histogram represents the distribution of 28 sBzKs 
with a spectroscopic redshift given in IDaddi et al.1 
(|2004l ). The solid curve indicates a Gaussian dis- 
tribution function with z c — 1.9 and a z = 0.35 
used to calculate the correlation length. The two 
Gaussians with a z = 0.25 and 0.45 are plotted by 
the dotted line to show the adopted uncertainty 
in the redshift distribution. The amplitude of the 
Gaussians has been normalized so that their areas 
be equal to the area of the histogram. 



the adopted uncertainty in a z (±0.1). 

Figure [Jj shows the correlation lengths for our 
sBzK sample and the six samples of K06. 1 It is 
found from Figure [7] that ro increases rapidly with 
increasing K brightness. The brightest sample of 
K06 has ro ~ 20 /i _1 Mpc, which is six times larger 
than that of ours. 



lr The errors in ro of the samples of K06 are not very large 
in spite of the small sample sizes, because brighter sBzKs 
are clustered more strongly and thus have stronger signals. 
This in turn means that larger samples are required to 
measure the angular correlation function for fainter objects 
with the same uncertainty. 
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Tabic 2: Numbers, ACF amplitudes, correlation lengths, and dark-halo masses of sBzKs in our sample and 
K06's 



Field N A[arcsec u «] r [h~ l Mpc] (M DH )[M S ] 



,+0.6 
-0.7 
91 +5.2 
zl -6.2 
io+4.5 
i8 -5.4 
14 +4.4 



-2.3 
+0.9 
-0.9 
1 1+O.8 

J-^-o.e 
6 4+ 46 

°-^-4.8 

1 q+0-8 

x - a -0.8 

^■ 8 -2.0 
5+1.4 
-1.2 



JU DH 



x - u -0.9 

8.oig:8 

4 2 +4 - 8 

+0.6 
0.6 
+ 1.7 
1.3 
+9.0 



[M 



OJ 



this work 
K06 



SDF 
Daddi-F 



Deep3A-F 



23.2 
20.4 
20.7 
21.1 
20.7 
21.4 
21.9 



1092 

21 

43 

92 

27 
129 
387 



0.58±0.13 
16.8±6.86 
13.2±5.26 
8.05±4.63 
20.6±4.06 
4.69±2.20 
3.46±1.18 



3.2] 



23 
10 
8.5 



5.7 
+4.1 
4.7 
+2.8 
3.4 
+2.0 
-2.3 



2. 

1.6 i 



1.8 



10 11 
10 14 
10 14 
10 13 
10 14 
10 13 
10 13 



1.5 
1.6 
9.0 



10 11 
10 14 
10 13 
10 13 
10 14 
10 13 
10 12 



5. Discussion 

5.1. Masses of Dark Haloes Hosting sBzKs 

The standard CDM model predicts that at any 
redshifts more massi ve dark haloes are o n average 
more clustered (e.g-. lMo fc White! [20021) . We use 
the clustering strengths of our sBzK sample and 
those of K06 to infer the masses of their hosting 
haloes in the following manner. First, for an sBzK 
sample with a given limiting magnitude, we derive 
the amplitude of the spatial correlation function at 
8 h^ 1 Mpc, ^ S BzK(8/i _1 Mpc), from the observed 
ro assuming £ S BzK oc r~ 18 . We then calculate the 
effective bias parameter, & e ff, from 



b c s = 



3 zK (8fe~ 1 Mpc) 



(1) 



where £dm is the predicted correlation function of 
dark matter at z = 2. Finally, we use the analytic 



formu la of dark-halo biasing given in ISheth et al 



(2001) to obtain the mass of dark haloes hosting 



the sBzKs. 

We consider two definitions of halo mass, (i) 
average mass and (ii) minimum mass. Average 
mass, (M DH ), is defined by: 



b eS = &«M DH », 



(2) 



where 6(Mdh) is the bias parameter of haloes with 
mass Mdh- We regard (Mdh) as the typical halo 
mass of the sample. Minimum mass, M™^ 11 , is de- 
fined through: 



DH 



f™ min n(M DH )dAf D H 



(3) 



where n(AfDH)<iMDH is the mass function of dark 
haloes. In this definition, M^ 1 corresponds to 
the mass of haloes hosting the faintest galaxies 
in the sample (In this paper we assume that bias 
increases monotonically with galaxy brightness). 

Our definition of M™^ can be regarded as 
a simplification of the halo model approach 
which uses three paramet ers including M"™ 
to model galaxy clustering (IBerlind fc W einberg 



2002: Bullock et al. 2002: Moustakas fc Somerville 



2002: lHamana et al.ll2004l ): the other two param- 
eters are used to model the mass dependence of 
the number of galaxies per halo. In this study, we 
do not adopt this approach, but place a constraint 
only on Afg 1 ^ 11 for simplicity since this approach 
requires an accurately measured ACF shape with 
small angular bin sizes. 

The average dark- halo mass of our sBzK sample 
is estimated to be (Mdh) — 3 x 10 11 M©. In con- 
trast, the sBzKs of K06 samples are found to be 
hosted by massive haloes with (A/dh) — 2 x 10 13 
- 2 x 10 14 M Q . 

The difference between (Mdh) and M™^ 11 is less 
than factor 3. This is because fainter objects are 
more numerous, thus contributing more to the 
overall clustering of the sample. For the readers' 
reference, The dotted horizontal lines of Figure [7] 
indicate M^ii as a function of r . 

These results show that sBzKs reside in dark 
haloes with a mass range of as wide as three or- 
ders of magnitude, 10 11_14 M Q , with fainter ob- 
jects being found in less massive haloes. Typical 
haloes hosting faint sBzKs (K < 23.2) are several 
times less massive than that of the Milky Way 
(M DH ~ (1 - 2) x 10 12 M Q : e.g.. ISakamoto et al 
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Fig. 7. — Correlation length, ro, as a function of 
limiting K total magnitude. The filled circle shows 
the measurement for our sample, while the filled 
triangles correspond to the six samples with dif- 
ferent limiting magnitudes over K = 20.4 and 21.9 
given in K06. The dotted horizontal lines indicate 
the minimum mass of dark haloes inferred from 
To- 



(l2003l) 1 while those hosting bright sBzKs with lim- 
iting magnitudes of K = 20.4 to 21.9 have masses 
comparable to those of present-day galaxy groups 
and clusters. 

5.2. Ratio of Stellar Mass to Dark Halo 
Mass 



We estimate the minimum stellar mass, which 
corresponds to that of the faintest galaxies in the 
sample, for our sample and the six samples of K06 
from the K magnitude and z' — K color of the 
faintest object i n each sample , using equations (6) 
and (7) given in iDaddi et all (|2004l ). These equa- 
tions were derived using bright sBzKs which have 
accurate stellar masses base d on a fit of mode l 
spectra to multicolor data. IDaddi et al.l |2004) 



found the typical error in mass estimates from 
these equations to be AlogM* = 0.2. 

Figure [5] plots the minimum stellar mass of 
sBzKs against the minimum mass of their host- 



ing haloes. The minimum stellar mass of our 
sample is Mf ln ~ 1 x 10 10 M Q , while those of 
K06 are three to ten times higher. On the other 
hand, the minimum dark-halo mass of our sample 
is ~ 1 x 10 n A/ o , while those of K06 are two 

to three orders of magnitude higher. It is thus con- 
cluded that dark-halo mass increases more rapidly 
than the stellar mass of individual sBzKs residing 
in the halo. 

We note that the stellar-mass estimation for our 
sample is less reliab le than that for K 06. Equa- 
tions (6) and (7) in IDaddi et all (|2004l ) have not 
been tested for faint sBzKs like those in our sam- 
ple, which may have systematically different stel- 
lar populations and dust extinctions. The ob- 
served trend of increasing M™™/M™ in with M™ in 
will disappear if the M™ m of our sample is over- 
estimated by more than factor ten, although such 
a large overestimation seems to be unlikely. 




10» 10'° 10" 10" 10" 

minimum stellar mass [M ] 

Fig. 8. — Minimum dark-halo mass plotted 
against minimum stellar mass for sBzKs. The 
filled circle indicates our sample, while the filled 
triangles correspond to the six samples with differ- 
ent limiting magnitudes given in K06. Minimum 
stellar masses correspond to the limiting magni- 
tudes of the samples. 
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5.3. Halo Occupation Number 

The halo occupation number of a given class of 
galaxies is defined as the average number of this 
class of galaxies hosted in a single dark halo. We 
estimate the halo occupation number of sBzKs for 
two magnitude limits, K < 23.2 and K < 21.9, us- 
ing our sample and the faintest sample of K06, and 
then examine whether or not occupation number 
depends on luminosity. 

The halo occupation number for sBzKs with 
K = 23.2 is estimated by dividing the number 
density of our sBzKs by that of dark haloes host- 
ing them, i.e., haloes more massive than Afgjj 1 . 
We have 1092 objects in a survey volume of 7.1 x 
10 5 Mpc 3 . This leads to a number density of 
2.4 x 10~ 3 Mpc~ 3 after correction for detection 
completeness of a modest amount. The num- 
ber density of dark haloes is calculated to be 
1.7 x 10~ 2 Mpc~ 3 from the an alytic formula of th e 
halo mass function given by ISheth et al. ( 2001 ) . 
Thus, the halo occupation number of sBzKs with 
K < 23.2 is estimated to be ~ 0.1. 

Similarly, we use the faintest sample of K06 
with K < 21.9 to estimate the halo occupation 
number to be ~ 4. 2 Here, the number of 
sBzKs in this sample is 387, the survey volume 
is 3.2 x 10 _4 Mpc -3 , and the number density of 
the hosting haloes is 7.5 x 10~ 5 Mpc~ 3 . 

These two estimates show that the occupation 
number of sBzKs increases strongly with increas- 
ing halo mass. Haloes with Mdh ^ 10 13 M Q have 
several sBzKs on average, while only one tenth (on 
average) of haloes with Mdh ^ 10 n M Q host an 
sBzK. It is reasonable that haloes with > 10 13 M Q 
host multiple sBzKs, since they are as massive as 
present-day groups and clusters of galaxies. On 
the other hand, the very small occupation number 
for haloes with > lO n M may imply that they 
are not massive enough to always feed a galaxy 
more luminous than K — 23.2. 



2 We do not use the brighter five samples of K06, because 
sBzKs in them are hosted by more massive haloes, whose 
number densities, and thus occupation numbers, have in- 
creasingly large errors owing to a steep slope of the mass 
function. 



5.4. Comparison of ro between sBzKs and 
Other Objects 

Figure [5] compares the correlation lengths of 
sBzKs with those of various types of objects from 
present-day early-type galaxies and clusters of 
galaxies, up to galaxies at z = 4. The five curves 
show the correlation lengths of dark haloes with 
five fixed values as a function of redshift. 

Haloes with a fixed minimum mass have larger 
ro at higher z because b increases with z for any 
Mdh- 

It is found from Figure [9] that at z ~ 2 faint 
(K < 23.2) sBzKs and BX/BM galaxies have sim- 
ilar correlation lengths and thus similar dark-halo 
masses of the order of ~ 1O 11 M0. This result 
i s in a ccord with the recent claim by iReddv et al.l 
(120051) that a high fraction (70 - 80%) of BXs/BMs 
with 22.3 < K < 22.8 are also selected as sBzKs 
and that the fraction increases with decreasing K 
brightness. Therefore, faint sBzKs and BXs/BMs 
are largely overlapping subsets of star-forming 
galaxies at z ~ 2 hosted by relatively low-mass 
haloes. In Figured! LBGs at z = 3 - 4 have also 
similar masses to faint sBzKs, while haloes host- 
ing QSOs are on average an order of magnitude 
more massive. 

Brighter sBzKs with K ~ 21 - 22 are found to 
have as large correlation lengths, ro ~ 10 - 15 hr 1 
Mpc, as DRGs at z ~ 2 - 3 and EROs at z ~ 1 - 
2. They are located in haloes as massive as the or- 
der of Mdh ~ 10 13 A/©. Their correlation lengths 
partly overlap with those of present-day luminous 
early-type galaxies and clusters. Large fractions of 
DRGs and EROs are known to be passive galax- 
ies, while sBzKs are star-forming galaxies. This 
implies that haloes with similar masses can host 
galaxies with extremely different star-formation 
properties. Brightest sBzKs with K < 21 are most 
strongly clustered among the high-z populations 
plotted here. Their correlation lengths are compa- 
rable to or even larger than those of present-day 
rich clusters. 

5.5. Present-day Descendants of sBzKs 

We use the extended Press-Sch echter formal- 
ism ( Bond et al. 1991k [Bower 1991) to predict the 
mass of the present-day descendants of dark haloes 
hosting sBzKs at z ~ 2. We find that haloes at 
z = 2 with a mass equal to the (Mdh) of faint 
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Fig. 9. — Correlation lengths of various types 
of objects over < z < 4. The filled sym- 
bols represent sBzKs; circle: our sample, trian- 
gles: K06. The other symbols indicate the mea- 
surements for other objects taken from the lit- 
erature; open circles: LBGs a t z ~ 3 (23.5 < 
R < 25.5, 
/ < 24.5, ' 



Adelberger et all (l2005l) : 20.0 



h003) and 



, Foucaud et al . 
{%■ < 25. 26. 27.5. buchi et all (|2005l )). open pen- 
tagons: BX/BM galaxies with 23.5 < R < 
25.5 at z ~ 2 (I Adelberger et all l2005f) . stars: 
QSOs at z ~ 0.5 — 2 (18.25 < bj < 20.85, 



Croom et al.l (|2005l )). open triangles: Distant Red 
Galaxies (DRGs) a t z ~ 2 — 3 (K Vesa . < 21, 



Qua dri et alj (120061): K < 20. 7 iFoucaud et al 



(|2006l ): K < 23.5. lGrazian et alj (|20btj|) ), open in- 
verted triangles: Extremely Red Obj ects (EROs 



at z 



2 (Kx 



_ . ifiga I ■ 

(120031): jfygga < 2 2 iRoche et al.l (|2003f) : if Vo ga < 
18.4. [Brown et all (|2005h ). open sq uares: present 



< 24, Daddi et al 



day l uminous early-type galaxies (jOverzier et al 
20031 ). d ouble circles: pres ent-day clusters of 



galaxies (Ba hcall et al.l 120031 ). The five curves 



show the correlation length of dark haloes with 
five different minimum masses as labeled. 



{K < 23.2) sBzKs (2.8 x 10 n M Q ) will become 
haloes with masses of (3.7 - 10.0) x 10 n M Q at 
z = (68% range of the distribution function; the 
extended Press-Schechter formalism can predict 



descendant masses only in a statistical sense). 3 
Thus, the present-day descendants of the typ- 
ical haloes of faint sBzKs are likely to be less 
massive than that of the Milky Way. They may 
host sub-i* galaxies in nearby fields. On the ba- 
sis of a clustering ana lysis of BX/BM galaxies, 
lAdelberger et alJ (|2005l ) concluded that BXs/BMs 
are progenitors of normal ellipticals in the local 
universe. The correlation length they obtained is 
r = 4.2 ± 0.5 h^Mpc for BXs and r Q = 4.5 ± 0.6 
/i _1 Mpc for BMs. These values correspond to dark 
haloes several times more massive than those of 
our faint sBzKs. Thus, their conclusion does not 
seem to conflict with our result obtained here. 

On the other hand, haloes with a mass equal 
to the (M DH ) of K < 20 .4 sBzKs in K06 (1.6 x 
10 14 M o ) are predicted to have (4.3-9.3) x 1O 14 M 
at present. These masses are comparable to those 
of most massive clusters like Coma. The num- 
ber density of sBzKs with K < 20.4, (1 - 2) x 
10 -5 Mpc~ 3 , is close to that of the present-day 
rich clusters mor e massive than 2 x 1O 14 M , ~ 
1 x 10~ 5 Mpc- 3 (|Rines et al.l 120061 ). This sug- 
gests that brightest sBzKs are ancestors of central 
galaxies in rich clusters seen at present. 

These calculations suggest that sBzKs evolve 
into galaxies over a wide range of mass (or lu- 
minosity) in a variety of environment at z = 0, 
depending on their apparent K brightness. 

6. Conclusions 

In this paper, we have studied clustering prop- 
erties of star-forming BzK galaxies (sBzKs) at 
1-4 ^ z ^ 2.5 over a wide range of K brightness 
(K < 23.2). 

We have used deep multi-color data of 180 
arcmin 2 in the Subaru Deep Field to construct a 
sample of 1092 faint (K < 23.2) sBzK galaxies. 
We have derived the angular correlation function 
(ACF) of the sBzKs, and measured its amplitude 
to be A = (0.58 ± 0.13) [arcsec 8 ] by fitting a 
power law of w(8) oc 6~ Q - S to the data. We have 
then transformed the ACF amplitude into the cor- 
relation length assuming a Gaussian redshift dis- 
tribution of z c = 1.9 and a z = 0.35 ± 0.1, and ob- 
tained ro = 3.2^Q 7 h~ x Mpc. We have not been 



3 Dr. Takashi Hamana kindly provided the code to calculate 
the distribution function. 
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able to derive the ACF of passive BzKs in our 
data because of the small sample size. We infer 
from the correlation length that our sBzKs reside 
in haloes with an average mass of 2.8 x 10 11 Mq 
and a minimum mass of 1.0 x 10 11 Mq. 

We then have inferred dark-halo masses for six 
bright (K < 21.9) sBzK samples of K06 in the 
same manner. Combining our data with those of 
K06, we have examined how the mass of host- 
ing dark haloes depends on the K luminosity of 
sBzKs. We have found that the mass of dark 
haloes rapidly increases with K brightness of indi- 
vidual sBzKs; K06's sBzKs are brighter than ours 
by only up to ~ 10 times, but they are hosted 
by haloes two to three orders of magnitude more 
massive than the haloes of our sBzKs. The halo 
occupation number, the number of BzKs hosted 
in a dark halo, is found to be higher for brighter 
sBzKs. 

The correlation lengths of sBzKs have been 
compared with those of various types of objects 
up to z = 4. We have found that faint (K < 23.2) 
sBzKs have similar r$ values to BX/BM galaxies, 
which are optically-selected star-forming galaxies 
at z <~ 2. We argue that these two types of galaxies 
are similar populations hosted by relatively low- 
mass haloes with the order of Mdh ~ 10 11 M Q . 
On the other hand, sBzKs with K ~ 21 - 22, 
EROs, and DRGs arc found to reside in massive 
haloes with the order of ~ 10 13 M Q . The clustering 
of brightest sBzKs with K < 21 is the strongest 
among the high-z populations. 

Finally, we have predicted present-day descen- 
dants of haloes hosting sBzKs using the extended 
Press-Schechter formalism. Descendants are found 
to span a wide range of mass, depending on the 
K brightness of sBzKs in them. Typical descen- 
dants of haloes hosting sBzKs with K < 23.2 are 
less massive than the Milky Way, while those for 
K < 20.4 sBzKs are as massive as richest clusters. 

We wish to thank Takashi Hamana for provid- 
ing his code to compute the mass growth of dark 
haloes using the extended Press-Schechter formal- 
ism. We are very grateful to Emanuele Daddi 
for kindly providing new data of the redshift dis- 
tribution of 81 sBzKs prior to publication. We 
also thank an anonymous referee for useful com- 
ments which have greatly improved the paper. 
The United Kingdom Infrared Telescope is oper- 



ated by the Joint Astronomy Centre on behalf of 
the U.K. Particle Physics and Astronomy Council. 
The WFCAM data were reduced on the general 
common-use computer system at the Astronomy 
Data Center, ADC, of the National Astronomical 
Observatory of Japan. 
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